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The anomalously strong scarring of wavefunctions found in 
numerical studies of quantum wells in a tilted magnetic field 
is shown to be due to special properties of the classical dy- 
namics of this system. A certain subset of periodic orbits are 
identified which are nearly stable over a very large interval of 
variation of the classical dynamics; only this subset are found 
to exhibit strong scarring. Semiclassical arguments shed fur- 
ther light on why these orbits dominate the experimentally 
observed tunneling spectra. 



The localization of certain quantum wavefunctions in 
real-space along unstable classical periodic orbits illus- 
trates how quantum mechanics can violate the ergodic 
behavior expected from classical mechanics. Such wave- 
functions are conventionally termed "scars" jjj and their 
properties have been extensively studied by theorists of 
quantum chaos during the past decade [|]-[|. Recently 
an experimental system has been discovered and stud- 
ied HI! m which such scarred wavefunctions control to a 
large extent an observable physical property, the tunnel- 
ing current through a double-barrier GaAs-AlGaAs het- 
erostructure ("quantum well") under high bias. When a 
magnetic field is applied at an angle 9 with respect to 
the normal to the barriers (the electric field direction), 
the resulting dynamics makes a transition to chaos |^-|7|] 
as 9 is increased from zero. Calculations by Fromhold 
et al. Q on the system found many more scars than 
in any previously-studied quantum-chaotic hamiltonian, 
and that these scarred wavefunctions carried most of the 
the tunnel current when the system was resonant. In 
the initial experiments |5|,[| the level-broadening (due to 
optic phonon emission) was too large to observe the res- 
onances due to individual levels. However in a later ex- 
periment ||] this was done, albeit at such low quantum 
numbers that the concept of scarring becomes less mean- 
ingful Jll[] . From extensive numerical J|,|],|l2],[l5) work we 
know: 1) Quantum states scarred by the same periodic 
orbit appear over a wide range of variation of the classical 
dynamics, in contrast to typical systems (e.g. billiards). 
2) The scars arise from only a few of the many unstable 
short periodic orbits in the tilted well. 3) These scars 
carry most of the resonant tunneling current for 9 > 15°. 
In this Letter we will present a theory to explain why 
only certain orbits scar the wavefunctions and why these 
scars persist as the classical dynamics changes substan- 
tially. The theory also sheds further light on why these 



scarred states dominate the tunneling current. 

First, we recall why in typical chaotic systems scarred 
states are relatively rare. Generally scarred states are 
associated with orbits which are not too unstable; a suf- 
ficient condition for strong scarring is that XT < 1, where 
A is the largest instability (Lyapunov) exponent associ- 
ated with the orbit. Typically in chaotic hamiltonian 
systems, unstable periodic orbits appear as marginally 
stable orbits at bifurcations and become monotonically 
more unstable as the classical parameter (e.g. energy) 
driving the system to chaos is increased. Therefore such 
orbits only scar over the small interval of classical param- 
eter space when they are close to stability. In a recent 
detailed study of the classical dynamics of the tilted well 
|l5| , we have shown that in this system a subset of the 
short periodic orbits behave completely differently. They 
exist only for a finite interval as the classical parameters 
(energy e, magnetic field, B and electric field E) are var- 
ied and are "pinned" near marginal stability for their 
entire interval of existence. Therefore these orbits are re- 
sponsible for the strong scarring seen, and indeed all the 
scars found numerically correspond to this subset of the 
short orbits. We now present this argument in detail. 

We model the system by two infinite potential barriers 
corresponding to the x — y planes at z = —d (the emitter) 
and z = (the collector), with an electric field E = Ez, 
and a "tilted" magnetic field B = B cos 9z + B sin 9y 
for — d < z < 0. The periodic orbit theory for this 
model is rather involved ]lq ] and we only sketch the 
most salient features here. First, the classical hamil- 
tonian can be rescaled jl5| so that the dynamics only 
depends on two dimensionless parameters: (3 = 2Bvq/E 
and 7 = e/eV where V is the voltage across the well 
and vo = (2e/m*) 1 ' 2 is the velocity corresponding to 
the total injection energy. In the experiments 7 m 1.17 
is constant to a good approximation |Te|| and so /3, the 
scaled magnetic field, is the single relevant variable. Our 
semiclassical analysis |1| expresses the tunneling current 
in terms of periodic orbits in the well which connect the 
emitter and collector barriers (we term these "emitter" 
orbits, those which don't reach the emitter "collector" 
orbits). We focus on the experimentally relevant situa- 
tion in which 7 is only slightly greater than unity, so that 
many of the periodic orbits are collector orbits and may 
be neglected. Optic phonon emission produces a tem- 
poral cut-off which allows only short periodic orbits to 
produce structure in the tunneling spectra. 

Emitter orbits have the following properties [ fT5| |, An 
orbit which collides with the collector n times (period-n 
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in the collector Poincare map) can collide with the emit- 
ter m times, where to = 1,2, ... ,n. Hence we denote 
the orbits as (m, n). All emitter orbits (except for a sin- 
gle (1, 1) "traversing" orbit) only exist above a threshold 
value of (3 = (3 C \ and cease to exist at a higher value 
of [3 = [3 C 2 by an inverse bifurcation (usually a tangent 
bifurcation (TB)). While the "death" of the orbits fol- 
low the generic rules of hamiltonian bifurcation theory 
Jbij , all relevant orbits are "born" in a new kind of bi- 
furcation, which we refer to as a cusp bifurcation (CB). 
CB's have non-generic properties since they appear on a 
closed curve in the surface of section (SOS), where the 
the Poincare map describing the dynamics is nonana- 
lytic. This "critical" curve separates initial conditions 
at the collector barrier which will reach the emitter on 
the next try from those which will not. Hence orbits 
originating just within the curve will receive a "kick" at 
the emitter, while those just outside will not. This leads 
to a discontinuity in the stability matrix M of any peri- 
odic orbit coresponding to a fixed point which crosses the 
boundary. In particular, all CB's occur by the simultane- 
ous appearance of two new orbits, which infinitesimally 
above threshold differ by one point of contact with the 
emitter (this may correspond to either one or two fewer 
collisions). We have shown jl5| that the orbit which 
reaches the emitter more times has diverging stability 
{\TrM\ — > oo ) at (3 c i while the other orbit in the pair can 
be either stable (\TrM\ < 2) or unstable (\TrM\ > 2) at 
Pd. The latter case, in which two unstable orbits appear 
simultaneously is forbidden for generic hamiltonian sys- 
tems 14 1 . Continuity arguments and numerical results 
[ p5| imply that the partner in a cusp bifurcation with 
fewer emitter collisions will have |TrM| ~ 2, i.e. will 
be born near marginal stability. Hence such orbits, be- 
ing born near marginal stability and being required to 
return to marginal stability when they disappear in an 
inverse tangent bifurcation, remain only slightly unstable 
for their entire interval of existence, unlike unstable peri- 
odic orbits of typical chaotic systems. It is this subset of 
the orbits which scar strongly the quantum states over a 
large variation of the classical parameter (3. 

Bifurcation and stability diagrams illustrating this be- 
havior are shown in Fig. 1 for the case of period-two 
and period-three orbits. Among the period-2 orbits the 
orbit denoted (1, 2)2 fits our criteria. At 9 = 28° it is 
born in a CB with the higher connectivity orbit (2, 2) + 
at j3 ~ 4.0 and dies in an inverse TB with the orbit 
(0, 2) _ at (3 ~ 7.0. This period-two orbit and another 
topologically similar orbit (not shown) which appears at 
slightly higher value of {3 account very well jl5| for the 
peak-doubling regions observed at 9 = 28° in the exper- 
iments of Muller et al. ||. This orbit was found to scar 
many wavefunctions in the work of Fromhold et al. @. 
It has a complicated evolution above 9 = 29°, involving 
an "exchange" bifurcation with the topologically simi- 
lar orbit just mentioned, the details of which are given 



elsewhere pjj|. Note however that in roughly the same 
(3 interval there are two other period- two emitter orbits, 
each with rapidly varying stability. Both born in CBs 
paired with an orbit with fewer collisions with the emit- 
ter, hence by our above reasoning are initially enormously 
unstable. Therefore they do not generate strong scars in 
the spectrum. 

A similar story holds for one of the eight period-three 
orbits which appear around (3 = 3.5 at 6 = 38°, the orbit 
we denote (1, 3)~ (Fig. 1). This orbit has been discussed 
previously |T^,|l6|,[l8| in connection with the observability 
of trifurcations in the data of ref . || . It is born in a CB as 
the partner of a (3, 3) orbit, remains near marginal stabil- 
ity for 3.2 < j3 < 4.4 and dies in a TB with the (1, 3)+. In 
the same interval there are several other unstable period- 
three emitter orbits which do not scar strongly. Finally, 
by the same reasoning we have found a (1, 5) orbit which 
scars strongly. Quantum states scarred by each of these 
orbits are shown in Fig. 2. 

Note that by our criteria the scarring orbits must al- 
ways be (to, n) orbits with to < n; e.g. (1,2) can scar 
strongly whereas (2,2) should not. On the other hand, it 
is easily shown [Jl5[ that as 9 — > the only emitter orbits 
are of the type (n, n). Therefore the interval of existence 
of the scarring orbits is small for small 9. Thus, for ex- 
ample, the period-three scarred states are unimportant 
for 9 < 20°. 

We have tested this argument quantitatively by analyz- 
ing the quantum states of the tilted well for scars of the 
three orbits shown in Fig. 2. By generating many spectra 
at different values of B, E we can search in the experi- 
mentally appropriate intervals of [3 with 1.1 < 7 < 1.2 
and systematically detect these scars. In Fig. 3 we plot 
a measure of scar stren gth versus action of the scarring 
orbit. As noted before p9|,§], the energies e n of scarred 
states satisfy an approximate Bohr-Sommerfeld quanti- 
zation rule, S(e n ) — (n + <fi)2Tth, so we expect and find a 
strong periodic modulation which turns on at the tangent 
bifurcation at which the orbit appears [po[ . 

Finally, we comment on the fact that these scarred 
states tend to dominate the tunneling current at large 
tilt angles. All of the orbits studied here and elsewhere 
HHjlJ] which scar strongly have only a single collision 
with the emitter barrier. Since both emitter and collector 
surfaces of section must be symmetric when v y — > —v v , 
such orbits must have v y = at the emitter barrier p2[ . 
Thus these orbits have unusually low transverse velocity 
at the emitter barrier compared to other orbits with the 
same periodicity in the collector Poincare map. Since 
the emitter wavefunction is primarily a superposition of 
the first few Landau levels, the source of tunneling cur- 
rent has low transverse velocity and couples very well to 
these scarred states. Therefore it is due to this specific 
feature of the emitter state that these scarred well states 
dominate the tunneling current; if the emitter state had 
large transverse momentum these states would be anti- 
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selected. Since scarred states are localized in real and 
phase space, they can lead to quasi-selection rules for 
tunneling, but only if they are localized in the correct 
regions to couple well to the input state. 
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FIG. 1. Bifurcation diagrams for the relevant period-two (a) and period-three (c) orbits. The horizontal axis is the scaled 
magnetic field, /3, the vertical axis is the v x coordinate of the fixed point(s) in the collector SOS indicated schematically in 
the insets, (b) and (d) are plots of the trace of the monodromy (stability) matrix for the corresponding orbits; Shaded area 
(|Tr[M] | < 2) denotes stability region. Orbits denoted (1, 2)2, (1, 3)~ (see Fig. 2) remain slightly unstable over a large variation 
of (3, leading to strong scarring. 
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FIG. 2. Examples of wavefunctions scarred by the unstable (1,2)2 orbit (a), (1,3)~ orbit (b) and a (1,5) orbit (c); y-z 
projections of orbits are superimposed. 
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FIG. 3. "Scar strength" ( Husimi function at the location of the fixed point of the periodic orbit at the emitter barrier, 
calculated for the normal derivative of the wavefunction) vs. the scaled action S(e„)/h of the corresponding orbit (e n is the 
energy of the corresponding eigenstate) for : (a) (1,2)2 orbit, (b) (1,3)~ orbit, (c) (1,5) orbit. The arrows indicate the values 
of f3 for the tangent bifurcation, which give birth to the periodic orbits, the peaks of the scar strength below these values 
are due to the "ghost effect" plpo) ]. Scaled actions below the bifurcation points were obtained by linear extrapolation of the 
(approximately linear) function S(e)/h 
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